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Abstract 

We studied theoretically the optical properties of Cqo fullerene microcrystals as a function of 
hydrostatic pressure with first-principles many-body theories. Calculations of the electronic prop- 
erties were done in the GW approximation. We computed electronic excited states in the crystal 
by diagonalizing the Bethe-Salpeter equation (BSE). Our results confirmed the existence of bound 
excitons in the crystal. Both the electronic gap and optical gap decrease continuously and non- 
linearly as pressure of up to 6 GPa is applied. As a result, the absorption spectrum shows strong 
redshift. We also obtained that "negative" pressure shows the opposite behavior: the gaps increase 
and the optical spectrum shifts toward the blue end of the spectrum. Negative pressure can be 
realized by adding cubane (CgHg) or other molecules with similar size to the interstitials of the mi- 
crocrystal. For the moderate lattice distortions studied here, we found that the optical properties 
of fullerene microcrystals with intercalated cubane are similar to the ones of an expanded undoped 
microcrystal. Based on these findings, we propose doped Cqq as active element in piezo-optical 
devices. 

PACS numbers: 1.48.-c,2.50.-p,71.35.Cc 
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I. INTRODUCTION 



Since its discovery, Ceo^ has been characterized as the most stable member in the series 
of fuUerenes, which are pure-carbon molecules with the shape of spherical shells^. Cgo 
can be produced economically and in abundance'^. Its chemical bonds have strong sp^ 
character, making the shell very stiff and at the same time free of dangling bonds. Together 
with clathrates and other carbon-based materials, fuUerenes are being actively investigated 
as building blocks for novel materials with unusual mechanical properties^!^. FuUerenes 
has remarkable properties, for example: Cqo doped with alkali and alkali-earth atoms is 
superconductor-i^; it has been claimed that Ceo can be heavily hydrogenated with up to 36 
hydrogen atoms per molecule, making it a promising material for hydrogen storage^. 

In pure form. Ceo crystallizes in a molecular solid (fuUerite), bound by weak forces between 
molecules. Its phase diagram is very rich, with a low-temperature face-centered cubic (FCC) 
phase, amorphous phases at intermediate temperature, and a diamond-like phase at high 
temperature. The FCC phase has rotational disorder and it is stable under pressure in excess 
of 15 GPa at room temperature^. Its energy gap is in the visible range, around 2 eV-^i^. The 
softness and stability of fuUerite could make it a good candidate for piezo-optical devices, 
where external pressure is applied reversibly and it modifies the optical response of the 
material. In addition, when fuUerite is heavily doped with molecules of appropriate size, it 
behaves as if it is under "negative pressure" . To that end, it is important to characterize 
the pressure dependence of the optical properties of this material. Extensive experimental 
work has been done on this direction^^ii^ii^'i^ii^'i^ii^. A review of the literature can be 
found in Reference joj. Theoretical analyses are not so extensive, mostly concentrated on 
characterization at zero pressure^ii^ii^S. In order to fill this vacuum, we present a systematic 
analysis of the optical properties of fuUerite at pressures ranging from zero to 6 GPa. We 
also investigate the optical response of the crystal at "negative" pressure, which can be 
realized in laboratory, for instance, by doping fuUerite with weakly interacting molecules 
such as cubane (CsHg). Numerical accuracy is essential, which is why we use state-of-the- 
art theories, namely many-body Green's function theories based on the GW-BSE approach. 

This article is organized as follows: we outline the theoretical framework in section [Tll 
That section is followed by a discussion of results at zero pressure in section llllt results 
at finite hydrostatic pressure in section IIVI and finally a description of the cubane-fuUerene 
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compound, which resembles fullerite with "negative pressure" , in section |Vl Finally, we 
conclude with some perspectives of future applications and a summary. 

II. THEORY 

The underlying electronic structure of fuUerene is determined using density-functional 
theory (DFT)^^. We use a plane- wave basis to solve the Kohn-Sham equations with cut-off 
in the kinetic energy of 50 Ry. Interactions involving valence electrons and core electrons are 
taken into account using norm- conserving pseudopotentials of the TrouUier-Martins type^i. 
We use the Perdew-Burke-Ernzerhof functional^^ (PBE) for exchange and correlation, based 
on the generalized-gradient approximation (GGA). It is well-known that DFT in the local- 
density approximation (LDA) or the PBE severely underestimates electronic gaps in general, 
making it unsuitable for detailed studies of optical properties of electronic systems. Accurate 
bandwidths and electronic energy gaps are calculated in a many-body framework within 
the GW approximation^^. In that approximation, the electron self-energy is computed by 
summing up Feynman diagrams to lowest order in the screened Coulomb interaction. At 
lowest order, the self-energy becomes a product between the one-electron Green's function 
G and the screened Coulomb interaction W , hence the name. We ignore vertex diagrams 
and we assume that Kohn-Sham eigenvalues and eigenvectors give a good approximation to 
the Green's function. Formally, the self-energy in space-energy representation is written as 

S(r, r'; E) = ^ f dE'e-^°^^'Go(r, r'; E - E')Wo{v, r'; E') , (1) 

where Go denotes the DFT-PBE Green's function and is the screened Coulomb interac- 
tion, related to the random phase approximation (RPA) dielectric function by: 

Wo{r, r'; E) = j dv"e-\v, r"; • (2) 

In the present formulation, the dielectric function is expanded in a basis of plane waves 
with cut-off 9.5 Ry. Its energy dependence is described by a generalized plasmon pole 
model^S. After the self-energy is computed, we diagonalize the quasi-particle Hamiltonian, 
H = HpBE + S — KrcT^'^i. Eigenvalues of that Hamiltonian provide the electronic band 
structure of the real material. This formulation is one of the simplest ab initio formula- 
tions of the GW approximation. Extensive applications of this formulation to a wide class 
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of carbon-based materials have shown it to predict electronic band gaps with an accuracy 
of 0.1 to 0.2 eV— In the specific case of fullerite, the first calculation of electronic gap 
within the GW approximation was consistent with direct /inverse photoemission spectra^. 
Owing to the fact that hydrostatic pressure on fullerite microcrystals does not affect their 
electronic properties besides an increase in intermolecular interactions, we expect our the- 
oretical methodology to be equally reliable in describing the pressure dependence of the 
electronic gap. Technical details about the theory can be found in review articles^^i^. 

The electronic band structure often does not give access to optical spectra because, after 
electron-hole pairs are excited, they interact and produce bound excitons, with energy lower 
than the electronic 

gap2^>26. 

We describe the dynamics of excitons by diagonalizing the 
Bethe-Salpeter equation (BSE) for electrons and holes. The BSE is an equation for the two- 
particle Green's function. Written as an eigenvalue equation, its solution gives the energy 
of optical excitations in the material: 

{E, - E,)Al, + Y: K'c'K'c' = ^'Al, . (3) 

v'c' 

where is the excitation energy of optical modes, indexed by s, and A^^ are the corre- 
sponding eigenvectors. K^P^, is the electron-hole interaction kernel, written in the basis of 
pair transitions. 

In the absence of electron-hole interactions {K = 0), each excitation energy is simply 
the difference between quasi-particle energies of electrons (Ec) and holes (Ey). The kernel 
K adds two types of interactions: an electrostatic interaction mediated again by Wo] and 
a repulsive exchange interaction between electron and hole, which is related to the fact 
that they can annihilate each other—. We follow the standard procedure to build and solve 
the BSE. We ignore the energy dependence of the interaction kernel and we assume the 
Tamm-Dancoff approximation^^- when computing K. Both approximations have been used 
extensively and they were shown to simplify considerably the numerical complexity, with 
little impact on numerical accuracy. This methodology has been presented in great detail 
elsewhere^i2^i2^. 

Although numerically expensive, the GW-BSE theory has been remarkably successful 
in predicting electronic and optical properties of real materials without the need for phe- 
nomenological parameters2ii2^. AH the approximations involved, such as the plasmon pole 
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model and the non-self-consistent assumption, are unambiguously defined. In addition, sp- 
bonded systems such as carbon-based nanostructures seem to be the ideal materials for this 
theory, owing to the fact that they have very weak correlation effects^^i^Ei^. 

III. FULLERITE AT ZERO PRESSURE 

The phase diagram of Ceo is extremely complex. At zero pressure and temperature, it 
crystallizes in a structure where the orientation of molecules is random but the molecules 
form a face-centered cubic (FCC) structure with lattice parameter around 14.2 A^. At 
room temperature and under mechanical pressure of 8 GPa, the crystals were found to 
polymerize in several phases, with substantial distortion of the cage^ii^. Since we are 
primarily concerned with hydrostatic pressure, we do not consider anisotropic pressure in this 
article. With increasing hydrostatic pressure, the lattice parameter decreases continuously 
according to Vinet equation of state^'^. The energy threshold of optical transmission also 
decreases^, following the reduction in lattice parameter. Other absorption edges are also 
known to redshift with applied pressure^. 

In our calculations, we apply pressure indirectly by fixing the lattice parameter and using 
Vinet equation to map lattice parameter into hydrostatic pressure^^: 



using a bulk modulus = 18.1 ± 1.8 GPa and its pressure derivative = 5.7 ± 0.6^^. The 
parameter x is the ratio between lattice parameters, x = a/ao. The lattice is built in the 
Pa3 (= cP12) structure, with one molecule per periodic cell. 

Figure [1] shows the calculated electronic and optical gaps for several choices of lattice 
parameter. At zero pressure, we obtain an electronic gap of 2.1 eV, in full agreement with 
previous work and compatible with photoemission and inverse photoemission datap^. Our 
DFT-PBE gap is 1.2 eV. The minimum gap is direct, around the crystallographic X point. 

The Ceo molecule has icosahedral symmetry, belonging to the Ih point group^^. Owing 
to its high symmetry, most molecular orbitals cluster in degenerate multiplets. The three 
highest occupied multiplets belong to symmetry representations denoted as II„, Gg, and 
(ordered from highest energy to lowest energy), with degeneracies 5, 4 and 5 respectively. 
The lowest unoccupied multiplet in molecular Ceo has symmetry Ti^, followed by a Ti^ 




(4) 
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multiplet, both with degeneracy 3. 

In fuUerite, each molecular multiplet originates a set of quasi-degenerate bands. The 
wavefunctions retain most of the shape of the molecular orbitals, so that they can still be 
labeled by symmetry representations of the molecular orbitals. The bandwidth of the 
quintuplet, at the top of the valence bands, is 0.5 eV. The next multiplet is a Tiu triplet, 
with approximately the same bandwidth. Within the GW theory, these bandwidths are 
slightly larger than the ones calculated with DFT-PBE. There are two major differences 
between band structures predicted with GW and DFT-PBE: (1) widening of the electronic 
gap, and (2) small stretch of bands according to the expressions: 



= E"p% X 1.2 + 0.6 eV 
E''^^^- = X 1.2 + 1.25 eV , (5) 

respectively for valence and conduction bands. In the equation above, the energies Eqw and 
EpBE are given with respect to the DFT-PBE valence band maximum. 

We determine the optical gap as the minimum excitation energy obtained after diago- 
nalizing the BSE. This gap at equilibrium lattice constant is calculated to be 1.7 eV. The 
oscillator strength associated to this excitation is very weak, owing to a molecular selection 
rule that prevents optical absorption from the Hu to the Ti^ multiplets. Significant absorp- 
tion is found around 2.2 eV, corresponding to Hu-Tig transitions, as shown on Figure [31 
The measured transmission edge is 1.9 eV— . This is compatible with our calculated results, 
considering the difficulties in determining the onset of absorption experimentally and the 
orientational disorder in the lattice, which is not included in our calculations. One of the 
earliest measurements of absorption spectra of crystalline Geo identified peaks at 2.0 eV, 
2.7 eV and 3.5 eV. Our calculations show peaks at 2.2 eV and 3.6 eV. 

The inset of Figured] shows the maximum exciton binding energy, defined as the difference 
between electronic gap and optical gap. The binding energy is high: around 0.4 eV. As 
discussed above, it arises primarily from the Goulomb attraction between electrons and 
holes, which is large compared with other solids because of a weak dielectric screening. The 
first bound exciton has well-defined Frenkel character, which is compatible with the fact 
that its binding energy is close to the electron and hole bandwidths. 

Figure [2] depicts the probability distribution of the electron given that the hole is fixed 

6 



on the surface of the central molecule. There are sharp maximums of probability on the 
central molecule, with more diffuse features in the neighbor molecules. 

In order to quantify the exciton radius, we have computed the integrated electron-hole 
probability and listed it on the third column of Table [H For the first bound exciton, the 
probability of locating electron and hole on the same molecule is 62%. The probability of 
locating the electron on any of the nearest neighbors molecule relative to the hole site in 
substantially smaller (30%), decreasing then to 2% if the electron is on any second nearest 
neighbor. Excitons with lower binding energy (and higher excitation energy) have more 
pronounced charge-transfer character, with the probability at nearest neighbor higher than 
the probability at the hole site. 

In order to address the validity of our calculations, based on the Pa3 lattice, with respect 
to the real, glassy crystal, we repeated the zero-pressure calculations with five different ori- 
entations of the molecule. As a result, the electronic gap fluctuated from 2.0 to 2.2 eV. That 
establishes an uncertainty in the determination of energy gaps arising from orientational 
disorder of the molecules. We find that orientational disorder affects similarly the electronic 
and optical gaps. The exciton binding energy fluctuates by less than 0.1 eV upon rotation of 
the molecular unit. Fine features in the absorption spectrum and in the density of states are 
smoothed out by molecular disorder while the broader features (energy position and width 
of major peaks) are very robust. 

IV. FULLERITE AT HYDROSTATIC PRESSURE 

Figure [1] shows that the electronic and optical gaps decrease continuously as an hydro- 
static pressure of up to 6 GPa is applied on crystalline Cqq. The overall decrease in electronic 
gap is 0.9 eV with pressure ranging from zero to 6 GPa. To our knowledge, the electronic gap 
at high hydrostatic pressure has not been measured yet. The optical gap {i.e., the excitation 
energy of the first bound exciton) decreases by 0.7 eV in the same pressure range. Since that 
exciton is optically inactive, the best comparison of optical activity as a function of pressure 
should be done following the position of the first peak in the absorption spectrum, on Figure 
m The peak moves from 2.2 eV to 1.75 eV in the pressure range from zero to 6 GPa. This 
is compatible with the first determinations of transmission edge as a function of pressure^: 
the transmission edge decreases from 1.9 eV (zero pressure) to 1.5 eV (5 GPa). 
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Figure [T] also shows that the profile of energy gap versus lattice parameter is not linear. 
A suitable model for the dependence of the gap with respect to pressure should take into 
account the behavior of dielectric screening for different amounts of intermolecular spacing 
and hence different amounts of overlap between molecular orbitals at different molecules. 
Snoke and collaborators^^ have proposed a phenomenological model for the gap. 

Meletov and Dolganov^ have also found a decrease in the optical g function of 

pressure. In their experiment, microcrystals of fuUerite were placed inside a diamond anvil 
cell, with pressure of up to 2.5 GPa. Several phenomena were observed in that experiment: 

1. At zero pressure, a low-energy line and two well pronounced lines in the absorption 
spectrum were found, labeled A (at 2.0 eV), B (2.7 eV) and C (3.5 eV) respectively. 
Line A is very weak and it could originate from transitions Hu — ^ Ti„, which gain 
finite oscillator strength from mixing with higher transitions. That interpretation is 
supported by our calculations, which indicate an onset of the line at 1.7 eV and very 
small but non-vanishing oscillator strength. 

2. Lines B and C have similar strength. It was found experimentally that optical activity 
migrates from C to B as the microcrystals are compressed. That effect is found in 
Figure [3], where we see enhancement of the peak at 2.2 eV and reduction of the peak 
at 3.5 eV, while both peaks redshift from zero to 3.4 GPa. Since we also see mixing 
between transitions Hu Tig and Hg — > Ti^, the major components of peaks B and C 
respectively, our calculations confirm the assumption that migration of optical activity 
is caused by mixing between different optical transitions^*^. 

3. The energy dependence of the measured absorption spectrum was reported to be 
weakly dependent on pressure in the pressure range from zero to 2.5 GPa^°. Figure 
[3] confirms that observation. At the next pressure value (6.1 GPa), the two peaks 
merge into an asymmetric wide peak. That indicates that bands derived from different 
molecular multiplets start to overlap, as shown in Figure HI 

Hydrostatic pressure also modifies the character of bound excitons, making them more 
delocalized. Comparing the distribution of probabilities at zero pressure (lattice parameter 
14.2 A) and 3.4 GPa (lattice parameter 13.6 A), Table [T] shows that the first bound exciton 
becomes substantially more delocalized, with a correlation radius between the first and third 
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nearest-neighbor distances. We believe that two mechanisms contribute to delocahzation: 
apphed pressure increases the overlap of molecular orbitals on different molecules, thus 
increasing the probability of one electron moving from one molecule to its neighbor; and 
pressure also increases the mixing between bands, particularly between transition —>■ Ti^ 
and higher transitions. 

As in the zero-pressure regime, we use the Pa3 lattice to perform calculations at hy- 
drostatic pressure, with no orientational disorder. Our estimates of the impact of disorder 
on the energy gap, mentioned at the end of the previous section, also apply to the regime 
of finite pressure. Since our calculations do not contain accurate van der Waals forces, we 
have not attempted to investigate the emergence of different glassy phases as a function 
of pressure. Including accurate van der Waals forces would remove the inaccuracy of the 
calculated gaps with respect to orientational disorder in fuUerite. 

V. FULLERITE WITH INTERCALATED MOLECULES 

FuUerene Cgo is very stable, which favors the engineering of microcrystals with interca- 
lated molecules. At equilibrium, the FCC crystal has two large types of voids: an octahedral 
site with radius 3.5 A and a tetrahedral site with radius 1.15 A. Isolated atoms and small 
molecules can be easily placed in one of those voids. Doped Ceo has very interesting prop- 
erties, for instance K3C60 is superconductor at 18 K^. CasCgo is superconductor at 8.4 K-. 
Those compounds also show significant electron transfer from dopant atom to cage. Doping 
fuUerite with wide-gap molecules produce different phenomena. Depending on the concentra- 
tion and symmetry of the dopant, it can lower the symmetry of the host crystal and enhance 
the oscillator strength of otherwise dark optical transitions of fuUerite. Highly-symmetric 
dopants are expected to produce less distortions in the host. In particular, cubane (CsHg) 
has been proposed as an ideal intercalator-ii. It has perfect cubic symmetry. If placed in 
an octahedral void, it will preserve the cubic symmetry of the lattice. With doping, the 
crystal is forced to expand isotropically in order to accommodate the extra molecules but 
no additional structural distortion is necessary. In addition, solid cubane is bound by weak 
van der Waals forces'^^, which means that cubane is not likely to segregate into clusters. 
The ionization potential of molecular cubane is 8.6 eV'^^. Its electron affinity is negative, 
indicating that it has an energy gap in the ultraviolet range. Since the band edges of Ceo 
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are inside the ones of cubane, cubane can be used to mimic negative hydrostatic pressure in 
fuUerite without altering the optical properties of the hoslr^. 

We built a lattice of cubane-fullerene with maximum doping by filling all octahedral voids 
in the FCC lattice with cubane. The lattice parameter is taken as 14.8 A, following exper- 
imental determination of the rotor-stator phase^i. The band structure of this compound 
around its energy gap is very similar to undoped fullerite with the same lattice constant. 
Cubane derived bands are found no less than 4 eV away from the gap. The electronic gap of 
fullerite with intercalated cubane (CgHg-Cgo), obtained from our GW calculations is 2.6 eV, 
similar to the electronic gap obtained for fullerite with the same lattice parameter (2.7 eV). 
The difference of 0.1 eV is close to the numerical precision of our calculations. CgHg-Cgo 
and pristine fullerite also have similar optical gaps: 2.0 eV and 1.85 eV. Figure shows that 
the absorption spectra of CsHg-Ceo and pristine fullerite differ from each other only above 
3.5 eV. These findings confirm that the effect of adding cubane is, by and large expand 
the lattice of fullerene molecules. All other phenomena in its electronic structure are direct 
consequences of lattice expansion. 

Other intercalants can also produce lattice expansion. Atoms of noble gases are good 
candidates, owing to their low reactivity and high energy gap. One shortcoming is that they 
are smaller than cubane. While these molecules will not expand the lattice significantly, 
a small change in gap can be measured. Therefore, fullerite could be used as a sensor of 
inert molecules. Significant expansion could be obtained by overdoping fullerite with several 
atoms per interstitial site. Other candidates are small molecules such as methane (CH4), 
hydrogen (H2) or nitrogen (N2). 

VI. SUMMARY AND PERSPECTIVES 

The results presented above show that several properties of fullerite, particularly the 
threshold of its optical absorption, can be tuned by applying pressure. By applying hy- 
drostatic pressure of up to 6 GPa, easily obtained in diamond anvil cell devices, the first 
peak of optical absorption redshifts from 2.2 eV (yellow-green) to 1.8 eV (red), thus making 
microcrystals less transparent. Similar reduction in the g function of pressure has 

been reported in alkali-doped fulleriteiS. Applications of this phenomenon are plenty. One 
of them is in piezo-optical sensors: one can put clean microcrystals of fullerite in an envi- 
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ronment under unknown pressure and infer the pressure by measuring their transmittance 
or absorbance. This is particularly useful if the microcrystals are part of a microdevice, 
subject to pressure gradients and where usual pressure gauges cannot be used. 

One can expand the range of colors where fullerite gauges operate by doping microcrys- 
tals with weakly interacting molecules. It has been shown experimentally that saturating 
microcrystals with cubane increases its lattice parameter-. Our results show that the lattice 
expansion produces a blueshift of the first absorption peak from 2.2 eV to around 2.6 eV. 
Other dopants can produce larger lattice expansion and hence larger blueshifts, depending 
on their size and concentration. 

In summary, we have done first-principles calculations of electronic and optical properties 
of fullerite in order to characterize their pressure dependence. Comparison between avail- 
able experimental data and our calculations at equilibrium lattice parameter show that our 
methodology predicts gaps with an accuracy of 0.1 to 0.2 eV. The absorption edge shifts 
toward the red end of the spectrum as we apply hydrostatic pressure of up to 6 GPa. There 
is little distortion in the electronic structure of the material in the pressure range inves- 
tigated. We have also confirmed earlier hypotheses that cubane-intercalated fullerite has 
optical properties very similar to fullerite with an artificial lattice expansion. These findings 
show that pure fullerite or fullerite with inert dopants can be used as active element in 
piezo-optical sensors. 

We would like to thank discussions with E. Schwegler, T. Oguitsu and H. Whitley for 
discussions. Research sponsored by the Division of Materials Sciences and Engineering BES, 
U.S. DOE under contract with UT-Battelle, LLC. Computational support was provided by 
the National Energy Research Scientific Computing Center. 
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FIG. 1: (Color online) Electronic gap (squares) and optical gap (crosses) calculated for fullerite 
as functions of lattice parameter. The equivalent hydrostatic pressure was obtained using Vinet 
equation, Equation HI The electronic gap was calculated within the GW approximation. The 
optical gap shown is the energy of the first excitation energy obtained from the BSE. It is a lower 
bound to the experimental optical gap since the first excitation has very low oscillator strength (see 
text). The inset shows the maximum exciton binding energy, i.e. difference between the electronic 
and optical gaps. 
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FIG. 2: (Color online) Isocontour plot of the electron probability distribution function of the first 
bound exciton provided that the hole is at a position where the highest occupied molecular orbital 
has maximum amplitude. The plot corresponds to fullerite at zero pressure. Only molecules up to 
second neighbor from the hole site are depicted in the figure. The isocontour shown corresponds 
to a value 10% lower than the maximum value of the probability distribution. This exciton is 
composed primarily by transitions in the Ti„ multiplet. 
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FIG. 3: (Color online) Imaginary part of the dielectric function for several choices of lattice param- 
eter. Vertical bars on each panel indicate the calculated optical and electronic gaps. An artificial 
Gaussian broadening of 0.02 eV was added to all absorption spectra. Sharp features in the spec- 
trum are expected to fade away with inclusion of rotational disorder. Line A in the measured 
spectrun>i^ (see text) is very weak to be visible. Line B is the first absorption line, at around 
2.2 eV at zero pressure. Line C is the second absorption line, at around 3.6 eV at zero pressure. 
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FIG. 4: (Color online) Density of states in fullerite for several choices of lattice parameter, obtained 
within the GW theory. Energies are defined with respect to the valence band maximum. A 
Gaussian broadening of 0.1 eV was added to all density distributions. The five major features 
in the density of states (well separated at zero and "negative" pressure) correspond to different 
molecular multiplets, from lower to higher energy: Hg+Gg (superimposed) , Hu, Tiu, Tig. 
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FIG. 5: Imaginary part of the dielectric function of CgHg-Ceo (a) and pure Ceo (b). An artificial 
Gaussian broadening of 0.02 eV was added to all absorption spectra. 
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Location of Electron 


a = 13.61 


Probability 
a = 14.21 


a = 14.81 


Hole's molecule 


15 % 


62 % 


84% 


I''* Nearest Neighbor 


25 % 


30 % 


13 % 


2nd Nearest Neighbor 


10 % 


2 % 


< 1% 


S'"'^ Nearest Neighbor 


25 % 


3 % 


< 1 % 



TABLE I: Electron probability distribution of the first bound exciton calculated as a function of 
the distance to the hole. The probability was integrated over each shell of molecules around the 
molecule that contains the hole. The integration was performed over a Wigner-Seitz cell centered 
on each molecule. Three different lattice parameters are shown: 13.6 1 (3.4 GPa), 14.2 1 (zero 
pressure) and 14.8 1 (negative pressure -1.6 GPa). 
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